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A cloud of cold A'^ two-level atoms driven by a resonant laser beam shows cooperative effects both in the 
scattered radiation field and in the radiation pressure force acting on the cloud center-of-mass. The induced 
dipoles synchronize and the scattered light presents superradiant and/or subradiant features. We present a 
quantum description of the process in terms of a master equation for the atomic density matrix in the scalar, 
Bom-Markov approximations, reduced to the single-excitation limit. From a perturbative approach for weak 
incident field, we derive from the master equation the effective Hamiltonian, valid in the linear regime. We 
discuss the validity of the driven timed Dicke ansatz and of a partial wave expansion for different optical 
thicknesses and we give analytical expressions for the scattered intensity and the radiation pressure force on 
the center of mass. We also derive an expression for collective suppression of the atomic excitation and the 
scattered light by these correlated dipoles. 
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1 Introduction 

Cooperative scattering by large collections of resonant atoms has been studied extensively for many years 
adopting either a classical or a quantum description. Often a quantum formalism is more convenient 
to describe the atom-light interaction. For instance, superradiant emission can be obtained when many 
independent atoms interact resonantly with photons, as studied in the seminal work by Dicke |JJ . However, 
many features can well be described with classical models for the atoms with a polarizability and for the 
light field. Many features of Dicke superradiance can thus be explained using classical theory. As atoms 
appear to be excellent systems to study any possible deviation from classical many body features, it is of 
general interest to understand and to monitor cooperative effects, even at a classical level, in a cloud of 
atoms. New intriguing effects can arise, when fluctuations due to the coupling with the vacuum modes 
can no longer be described by a classical field approach and the atoms can also become entangled during 
the cooperative scattering. Such cooperative scattering with atomic ensembles can appear in a number of 
experimental situations in free space |2| or in cavities [3|. 

Collective spontaneous emission at single excitation level |[T| has also been studied in resonant nuclear 
scattering of synchrotron radiation and has recently received growing interest with the study of single 
photon superradiance from N two-level atoms prepared by the absorption of a single photon |'6'-'9l. It 
has been shown that the photon is spontaneously emitted in the same direction as the incident field with 
a cooperative decay rate proportional to N and inversely proportional to the size of the atomic cloud ||8) . 
In a series of theoretical |101 and experimental flT,T2| studies, the authors have recently addressed the 
question of the quasi-resonant interaction of light with clouds of cold atoms, bridging the gap from single 
atom behavior, with granularity effects due to the discrete nature of the atomic distribution, to a 'mean- 
field' regime, where a continuous density distribution is the relevant description and leads to cooperative 
effects. In these works, the authors describe the collective atomic response under continuous excitation 
of a low intensity light field using an effective Hamiltonian model valid in the linear regime. From it, the 
average radiation pressure force and the scattered intensity have been derived, observing the modification 
from the single-atom values due to cooperativity 1 13 1. 
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Fig. 1 (color online) Experimental configuration: a cloud of two-level atoms (a) is driven by an incident laser detuned 
by Ao from the atomic resonance Wa, with wave vector ko (b). 



In this paper, we present an extended model based on a master equation in the single-excitation ap- 
proximation. Such model provides a more complete description of the linear regime, useful to investigate 
the distinctions between classical and quantum features, even though a complete description of quantum 
correlations in our system requires to take into account a larger number of excitations [14. 15J . In particu- 
lar, we show that by a first-order perturbation expansion of the master equation we recover the previously 
adopted linear model from classical optics. We also derive an analytical expression showing the coopera- 
tive suppression of atomic excitation and the related fluorescence. This effect bears some common features 
with the dipole blockade studied in Rydberg systems 1 16 ■ TsJ . In our case the long range l/r dipole-dipole 
coupling is at the origin of this important suppression even for dilute clouds of cold atoms. We note that 
this suppression occurs in the linear optics regime and is thus not related to a photon blockade. 



2 The master equation approach 

2. 1 The exact master equation for driven atoms 

The dynamics of a system of atoms driven by an external laser beam and undergoing cooperative re- 
emission into vacuum modes (see Fig. [T]) can be described by a master equation approach f 19''201. Let us 
consider a system of two-level atoms with transition frequency uja, positions Vj and excited decay time 
r. Each atom is described by the spin half angular momentum algebra, with 5*1 = \giji^ei\,S\ = |ei)((7i|, 
^\ = |ei)(ei| — l.9j)(.9j| satisfying the commutation relations \S\, S'i] = ^i3S\ and \S\, S*^] = =F2^yS'j_. 
The interaction Hamiltonian is 



N 



^-lAoi+iko-rj 



N 



Oke 



(1) 



j=l k 



where fto = dEo/h is the Rabi frequency of the classical incident field with amplitude Eq and wave 
vector ko, d is the dipole matrix element, Aq = ujq — oja is the pump-atom detuning, Ok is the photon 
annihilation operator with wave number k and frequency ujk = ck, gt = d[ujk/ {'2fieQVph)Y^'^ ^^'^ ^ph 
photon volume. We have assumed the rotating wave approximation (RWA) in the first term accounting for 
the interaction with the external laser, but not in the second term: there, the coupling between atoms and 
vacuum field modes is described in the scalar light approximation, where near field and polarization effects 
are neglected, since we are considering dilute clouds, with N{X/R)^ << 1 (R the system size). 
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It has been shown in 1 19 1 that the spontaneous emission properties can be conveniently described by a 
reduced master equation for the atomic system in the Bom-Markov approximation, given by 



dp 
dt 



.^0 



j E E {^SlpS^ - SXSip - pS^Si } . 



(2) 



where 



r cos(fco|ri — rj 
2 fco|r,-r,| 



(3) 



sin(fco|r.-r,|) ^ (4) 

Koki - Til 

and T = (PkQ / (27r?ieo ) ■ The first term in the master equation ^ describes the interaction with the external 
laser The second term describes the dipole-dipole interactions and arises from the virtual photon exchange 
between pairs of atoms. It becomes especially important at small interatomic distances, it is responsible 
for the collective Lamb shift |21 22| and plays an important role in the subradiant emission |(9] [23][24| . 
Finally, the third term of Eq. Q describes the cooperative emission. The Markov approximation ignores 
retardation effects and takes the long time limit, i.e. t ^ I/wq and t ^ max{|ri — rj\/c}. Hence, it 
requires that the system size is not too large, such that one photon travels through the atomic cloud faster 
than the characteristic cooperative emission time. 

The RWA approximation adopted in the present model requires more subtle arguments. As discussed 
in p9| , when applied to the second term of the Hamiltonian of Eq. ([T]i, the RWA consists in ignoring 
anti-resonant terms like aj^S*:^ and Ok'S'i which correspond to simultaneous creation or annihilation of 
a photon and atomic excitation (i.e. virtual transition). These terms are responsible of the shift of the 
ground state which contributes to the terms in Eq. (j2|. However, the master equation (|2ji has been 
obtained in 1 19] from the complete Hamiltonian ([T]l including the anti -resonant term, and making the RWA 
only on the master equation, neglecting the rapidly oscillating terms like S':^p5!j.e^*"°*. Hence, Eq. (j2jl 
obtained by making the RWA on the master equation rather than on the Hamiltonian does include the shift 
of the ground state. These remarks make clear that RWA on the Hamiltonian is not the same as RWA on 
the master equation and that one should make RWA on the final equations of motion. Since the counter- 
rotating terms as S'+pS'+e'^*'^"* are not important because T <C wq, Eq. ^ provides an accurate description 
of the interaction, including both the contributions of the real and virtual transitions to the frequency shift. 



2.2 Single-excitation approximation 

The master equation (|2]i has been used to describe, in the absence of the driving laser, the superradiant and 
subradiant decays from excited atoms. For a system confined in a volume whose size is much smaller than 
the radiation wavelength, and neglecting the dipole-dipole interaction terms A^, Dicke |TJ introduced the 
angular momentum states \ J, M) = Sym{|e . . .e;g . . . g)}, with M = — J, . . . , J. If the system is initially 
prepared in a symmetric (superradiant) state (with e.g. J = N/2), without any coupling between this state 
and the antisymmetric (subradiant) states. Note however that the presence of dipole-dipole interactions 
may induce coupling with subradiant states p3) . In our case, the presence of the driving laser provides the 
atom excitation which subsequently scatters and/or decays cooperatively. If the driving field is sufficiently 
weak or largely detuned from the atomic resonance, we can assume that the atomic system is weakly 
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excited, with at most one atom out of N excited. If the atoms are organized in a symmetric state, the decay 
rate is NT. Recently, it has been shown that an extended system of size i? ^ A (A = 27r/fco is the light 
wavelength) containing a single excited atom among N and prepared in the timed Dicke state ||4|[6l[8][25| 



1 ^ 

\TD) = -— ^cxp(jko • Tj) \gi, . . . ,6^,. . .,gN), (5) 



decays with a rate F^v oc NT / (koR)'^ . Also, an external field JIq drives the system mainly in a steady state 
of the form | ,10„11J 



where bo = (3A^/27r) / dz p(0, 0, z) cx N/ (koR)'^ is the resonant optical thickness along the propagation 
direction z of the driving laser and p(r) is the atomic density. However, a small fraction of the initial 
ground state is still coupled to the subradiant states, as discussed in [26 J . 

The previous works studying cooperative scattering by weakly excited atoms were based on a effective 
Hamiltonian model and brought us to a description equivalent to that of N classical linear dipoles driven 
by the external field 1 10, 13 |. Here, we go beyond the linear optics approximation using a master equation 
approach, still restricted to a single excitation. This restriction would be limited to the case of weak 
driving field, but may have a larger validity, for instance for Rydberg's atoms where some kind of blockade 
is provided |27J . 

By projecting Eq. (j2|l on the ground state \G) = \gi, . . . ,gN) and on the single-excitation states 
\i) = \gi, . . . , Ci, . . . , (?Ar), neglecting the states containing more than one excitation and defining pi o = 

{i\p\G) exp{iAot), pg,g = (G|p|G) and p,j = {i\p\j), we obtain 



dpG,G 

dt 

dpi,G 
dt 



dt 



.^0 



k,l 



lik 



k,l 



(7) 



E +i^ik) Pk,G, 

k^i 

no 



(8) 



y'"""'/OG,j -e-'*'''-""^p,,G) -i^A,kPk,j +i^AkjPi 



k^i 



E (^^kPk.j + IkjPt.k) ■ 



(9) 



We note that these equations still conserve the probability: pg.g + X^i Pi,i = 1- 



2.3 Perturbative solution 

We now show that the effective Hamiltonian approach can be obtained by a perturbative expansion of p 
in terms of the Rabi frequency fig of the external field. Assuming p = p'^^^ + p*^^' + . . . , if initially the 
system is in the ground state, the zero-order term yields p^gg ~ ^ pfc ~ pf*] ~ whereas at the 
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first order, Eqs. (|7]l-(|9]l reduce to: 



dt 
dt 



E(i) 



k,l 



E 



lik 



2 +i^ik)Pk,'g. 



.(1) 



dp. 



(1) 



dt 



EA (1) / (1) I (1) 

'^kjPi,k ~o2^ [likPk.j + P^M^k3 



(10) 

(11) 

(12) 



Since Pij{0) — 0, Eqs. (|lo| and (12i imply that at all times Pjj (t) = Pg g(^) ^ ^- '^'^^ remaining Eq 
|TT[ ) can be written defining Pi ~ and using Eqs. (jiji and as 

r ^ exp(ifco|r,, - rj. 



I ° 2/^' 2 



. . iko\r,-rk\ ' ''' 



(13) 



Furthermore, we observe that, since pi,j — {i\p\G){G\p\j) + 'Yl,k{Ap\^) {^\p\i) ^ the first non zero order 
for the dipole correlations corresponds to the second order for the field: p^^j* — Pj- q/Cq ^ = PiP*. 



2.4 Effective Hamiltonian 



Eq. (13 I can be expressed in the form of a Schrodinger equation 



dt 



= i?cff|*), 



where Y^) = a\g) + X^ili AK) ™d the effective Hamiltonian is 



(14) 



(15) 



where 



cos(fco|ri - r^l) .sin(fco|ri - rjl) 



(16) 



Assuming that at low saturation the system is weakly excited, so that awl, the projection of Eq. ( 14 1 over 



the excited states \i) leads to Eq. ( 13 i. Is has been shown that this equation describes also the temporal 



evolution of N harmonic oscillators driven by the scalar electric field radiation Eq, as predicted by classical 
linear optics [9 ,28]. It is important to notice that, even if the collective atomic state \^) is entangled, the 
knowledge of only the probability amplitudes Pi obtained from the linear equations ( [T3| is not by itself 
sufficient to detect entanglement, due to its classical nature. Conversely, entanglement could be observable 
in the solution of the master equation, even restricted to a single-excitation, since it contains correlations 
between atoms, and in particular between the ground state and the excited states. 

Finally, we notice that the exact master equation (|2]) can be written in terms of the effective Hamiltonian 



15 1, using Eqs. ( 14i, ( 15 i and p = |*)(^'|, as 



dp 

'dt 



YT.i^,sips\. 



(17) 



The last term in Eq. ( 17 1 describes the refilling of the ground state due to spontaneous decay of the excited 



states, and is necessary to preserve the density operator trace equal to unity. 
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2.5 Timed Dicke state 

Let us show how the steady state Eq. (j6]l provides an approximated solution of the single-excitation master 
equation. First, we observe that the external field couples the ground state \G) to the timed Dicke state 
\TD) defined in Eq. (jsjl. The matrix elements involving the TD state are 

PTD,G = (ra|p|G)--^5]e-^''«--V,,G, (18) 

j 

PTD,TD = (ra|p|ri?> = ^^^e-^''-(--^---Vj,™- (19) 

j m 

Their temporal evolutions are obtained from Eqs. (|7]l-(j9]) as 



= -« ^ [PTD.G - PG,TD) + 2^2^1klPl,k, (20) 

k I 



dt V 2 

1 



«Ao - - ptd.g - « — ^ — yPG,G - Ptd,td) 



dpTD,TD .^/N^O , . 

'* — ?, — [Pg,td — Ptd,g) 



dt 2 



- 2^ E E E e-*'-("^--^) h^kPk, + lk,P^,k) . (22) 

i j k 

These equations show explicitly the coupling induced by the laser between \G) and \TD). However, the 
dipole-dipole interactions and the cooperative decay favor the coupling also to all the other "subradiant" 
states |s), with s = 1, . . . , iV — 1, completing the single-excitation Hilbert subspace. We make the assump- 
tion that all the matrix elements among these states \s) can be neglected, supposing that their occupation 
probability is small compared to the timed Dicke state. In practice, we assume 

- {i\TD){TD\p\G) +J2{^\s){s\p\G) « ^e^'^-'-'pTD.G, (23) 



/N 
and 

p,,, = {t\TD){TD\p\TD){TD\j) + • • • « j^e'''°<''--''^'> Ptd,td- (24) 



By substituting these expressions in Eq. (20l-(22i we obtain 

dpG,G .VN^o 



dt 2 

dpTD.G 



{ptd,g — Pg,td) + ^nPtd,td, (25) 



dt 



•A 



PTD,G-i — ^ — \Pg,g - Ptd.td) , (26) 



dpTD,TD .VNno . , 

— — * — ^ — (Pg,td — Ptd,g) — nPtd,td, {^O 
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Atom Number Atom Number 

Fig. 2 (color online) Left: Ratio of tlie excited state population computed numerically and in the timed Dicke ap- 
proximation P/ Ptd- Right: Standard deviation ap. The cloud is Gaussian with fcofjj = 15. The contour plots show 
the iso-optical thickness lines 6(Ao) = &o/(l + 4Ao/r^) where bo = SN/{koa-R,)^ . This color-coded plot helps 
visualize the timed Dicke approximation validity region: 6(Ao) < 1. 



where A^r = Aq — and (T3\ 



1 AT r°° 



(29) 



are the cooperative decay rate and the collective Lamb shift, respectively. In Eqs. (28i and (29i we 
introduced the structure function 



1 ^ 

5Ar(k) = -^e^(''-''«)-"-^-. (30) 

and the average is taken over the total solid angle of emission of a photon with wave vector k at an angle 
6 with ko, where |k| — ko- Finally, the integral over n in Eq. ( [29) l is evaluated as a principal part. 

Eqs. ([25]l-([27ll have the form of Optical Bloch equations for a two-level system with collective states 
\G) and \TD), interacting with a collective Rabi frequency V NHq, detuning Ajv and Unewidth T^. The 
steady-state solution is given by. 



PTD,G ^ II ( OA I I ' ^^^^ 

1 + Sc \ 2/\n + iL n J 

pWto = ^(^), (32) 



where 



2Nnl 
4A2 ' 



(33) 



is the collective saturation parameter. We recover the previous linear result of Ref. pO) in the limit Sc ^ 1, 
with a steady-state amplitude of the excited state 

13^ « ( ^ e^k-"-- = ^e'''o r._ (34) 
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We have numerically verified that the timed Dicke ansatz of Eqs. ( [23] l and ( [24] i, or equivalently Eq. p4] i 
in the linear regime, is a good approximation of the system when 5(Ao) < 1, where b{Ao) = 6o/[l + 
(2Ao/r)^] is the optical thickness. Figure |2j shows a contour plot of the ratio P/Ptd (P — j \ Pj \ ^ 
and Ptd = I/^tdP, where j3j and jiro have been obtained from the numerical solution of Eqs.inJJi and 
from ( 34 1 respectively) as a function of N and Aq /F, for a Gaussian spherical distribution with parameter 
koCTR = 15. For all values of 6(Ao) < 1, the excited state population is well described by the timed 
Dicke ansatz and yields a ratio P/Ptd ~ 1- This result is confirmed by the analysis of the standard 

deviation ap = -y/(|/3p) - |(/3)|V|(/3)| of /3 = iSC^"-" (see right Figj2]where the contour plot of ct^ is 
also showed), ap quantifies the deviation of the excitation field from the 'mean-field' timed-Dicke state: 
this deviation becomes significant when 6(Ao) is larger than unity. Figj2] (right) suggests that for large 



optical thickness the homogeneous approximation assumed in the timed Dicke ansatz (see eq.(34i) is no 
more valid and a better description is demanded. 

2.6 Beyond the timed Dicke state approximation 

To account for the non-uniformity of the excitation within the cloud, the field (3 can be decomposed as 
a sum of waves that describe its spatial fluctuations. This approach has been considered in Ref. p9| , 
where Eq.lpjj) has been solved analytically for a continuous distribution with a Gaussian spherical profile, 
neglecting the cosine part of the exponential kernel and so the associated Collective Lamb shift A^ . Al- 
though the use of a such truncated kernel is not allowed when the decay of the excitation is observed pO) , 
it may still provide a reasonable approximation for small density and low optical thickness &(Ao). The 
stationary excitation amplitude obtained in ref. [29J using a partial wave expansion is 

= E 9A'"i^rni\ y U{kor)Pn{cos0), (35) 

where j„ are the spherical Bessel functions, P„ the Legendre polynomials, 9 the angle with respect to the 
laser wave vector ko and A„ = A^(7r/2cr^)^/^/„^i/2(o'^) exp(—o'^/2), where ct — koan Siud a fi is therms 
width of the Gaussian distribution. Eq. ( [35| ) accounts for non-homogeneity of the excitation probability 
density |/3(r)p. In this context, the timed Dicke expression (34i appears as a 'mean-field' approximation. 



which can be recovered assuming r(H-A„) ~ Tn inEq. (35l. Going further beyond, an exact solution for 



the continuous -density limit of Eq. ( 13 i with the exponential kernel has been derived using the Mie theory, 
although more mathematically demanding |31|. In that case, the excitation of an atom was calculated 
properly including both the incident and the phase-shifted radiation field from the other atoms. 



3 Observables 

3.1 Force on center of mass 

Cooperative effects can be investigated by a direct detection of the scattered photons. However, this mea- 
surement can be in general difficult, since for an extended atomic system the emission is strongly forward 
directed and the detector can be saturated by the incident laser The scattered radiation field detected at 
distance r and time t is the sum of the single fields scattered by the N atoms of position rj. 



dkl ^'^^o(t-\r-r,\/c) 

2ieo 



Eir,t) ^^Y. ^ ;„ \ ; ' SLit). (36) 



In the far field limit, |r — | w r — (r • rj)/r and 
2ieQr 



.2 N 
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where k = /co(r/f) and luq — cfco, so that the average intensity is 



I{r,t) = eoc{E\r,t)E{r,t)) - 



-ik-(rj-r„ 



(38) 



Alternatively, it is relatively easier to detect the cooperative effects by considering the radiation pressure 
force exerted on the atoms. If the atoms are sufficiently cold, it is experimentally possible to measure the 
atomic motion after their exposition to the incident laser beam. The radiation pressure force acting on the 
jth-atom is Fj = —Vy-.H 



0-] 



Fe, where 1 10 13 1 



-iAgt+iko-r, 



5; 



"A 



_ ■p(sclf) 



9 2^ 



- *A:oVn)e''=°''-^ + /i.e.} , 



(39) 
(40) 



where r, 



and Tjrn 



Faj and F^j result from the recoil received 



upon absorption of a photon from the pump and from the emission of a photon into a direction k, re- 
spectively. The emission force Fgj acting on the jth-atom has two contributions: a self-force F^^'^'^'' — 

^^^J2\k\=ko kS'^.S'i due to its own photon emission, and a contribution accounting for the coupling 
between the jth-atom and all the other atoms. Note that the dipole-dipole interactions can occur via a 
coupling to common vacuum modes of radiation. The interference terms in the total scattered field can 
leave a fingerprint on the forces acting on the atoms inside the cloud. This force has a term decreasing as 
1/rjm and one decreasing as Their average values on the single-excitation atomic states are 



(Fi) 



h.c.} 



(41) 
(42) 



2 ^ {korj^y 

m— 1 ^ ^ 

Notice that the self-force average to zero since the emission is isotropic. The force on the center-of mass 



of the atomic cloud, (F) = i'^/N)J2j(Pj) is of particular interest. FromEq. (41iand(42i, its component 
along the z axis of incidence of the laser is 




(43) 



where ji (z) — sin(z) / z'^ ~ cos(z) /z is the first order spherical Bessel function and Zj 



{Zj Zmj/rj 



In the timed-Dicke limit, using the approximations (23i, (24i and neglecting saturation, Eq. (43 i be 



= hkoT 



"o 



4A^ 



^ N 



N{{l-cose)\SNik)f)^ 



(44) 



The radiation pressure force can be influenced by different effects. On one side, the finite extent of the 
atomic cloud can produce strong forward oriented scattering. The balance between the momentum of the 
incident and scattered photons and the atoms indicate that for forward emission, the net recoil imprinted 
onto the atoms is vanishing, resulting in a reduction of the radiation force. A different contribution to the 
reduction of the radiation force can be seen in the prefactor of Eq. (44 1, which would appear even in the 



case of isotropic scattering (i.e. when (Fe) = 0). The importance of this prefactor can be understood from 
the cooperative coupling of several atoms into the same vacuum mode. The number of available modes 
for large spherical clouds can be estimated by N„i ~ {koR)^ (where R is the cloud's size), resulting in a 
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Atom number Atom number Atom number 



Fig. 3 (color online) Ratio of the cooperative to independent radiation pressure force as a function of the atom number 
A'^ and detuning Aq. The cloud is Gaussian with koUH — 10. The contour plots show the iso-optical thickness lines 
6(Ao) — 6o/(l + 4Ao/r^) where 6o = 'iN / {koa b)'^ . The left figure shows the numerical results computed from the 
effective Hamiltonian Eq. \\5\ and Eq. t he central figure stands for the analytical formula Eq. ( |45[ > and the right 
figure correspond to the modal expansion l|46^. 



number of atoms per mode scaling as N/N^ ~ N/ {koR)^. This scaling is conveniently related to the on- 
resonant optical thickness of the atomic cloud bg. For a Gaussian density distribution with root mean square 

size an, bo = 3N/{koaR)\ N{\SN\^)e,4, = l + bo/l2and N{{l-cos9)\SN\^)e,4> = bo/24{koaR)^ 1^. 

It is convenient to compare the cooperative radiation pressure force to the force acting on a single 
independent atom i^ind = 7ifcorrig/(r^ + Aq). The ratio of the cooperative radiation pressure force for a 
Gaussian cloud to the single-atom force in the timed Dicke limit can be written as (neglecting the collective 
Lamb shift) 



^md 4A2+(l+f|)^r2 



2 , 



1 ^° 



24{koaR) 



(45) 



One can also use the partial wave expansion to account for the inhomogeneity of the field /3, in which case 
the force ratio reads I 



^ T^ + Al ^ / (27i + l)A„(l + A„) _ (2n + 2)A„A„+i[4A§ + T^jl + A„)(l + A„+i)] 
F.„d N £-^UAg + r2(l + A„)2 [4A2+r2(l + A„)2][4A2 + r2(l + A„+i)2] 

(46) 

Fig. |3] shows the ratio of the cooperative to independent radiation pressure force as a function of the atom 
number N and detuning Aq for a Gaussian cloud with root mean square size fco^'i? — 10- The left figure 



shows the numerical results computed from the effective Hamiltonian Eq. ( 15 i and Eq. (43 i, the central 



one describes the timed-Dicke formula Eq. (|45]l and the right picture stands for the partial wave equation 



(46 1. While all methods predict a significant reduction of the force ratio for small detuning and large atom 
number, the partial wave approach yields better agreement with the numerical approach for Ao < F and 
N > 1000, despite the fact that the present partial wave approach has been limited to a sine kernel, not 
fully accounting for virtual-transition induced phase shifts BO). 



3.2 Dicke subradiance 

Dicke subradiance is the counterpart of superradiant emission and corresponds to the partial trapping of 
light due to destructive interferences. In a subradiant state, the atomic dipoles are arranged such that 
the macroscopic polarization of the cloud is small reducing the emission rate of the system. Subradiant 
emission has been previously observed for two ions p2| and also for the emission of a cloud of N atoms in 
a free space into a single radiation mode f33l. In a recent paper p6| , we showed that the system presented 
above is ideal to observe for the first time long photon storage into metastable subradiant states for N 
atoms in free space. 



In section 2.5 we saw that the laser pumps the system from the ground state |G) into the timed Dicke 



state \TD). Then the dipole-dipole interaction terms couple the timed Dicke state to the different other 
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Mode energy (r) 



100200300 



Fig. 4 (color online) Energies and decay rates of the modes of the system obtained by computing the eigenvalues of 
the effective Hamiltonian i^eff for Gaussian clouds with = 2000 atoms and koaa — 10, 30, 100 (feo ~ 60, 6.7, 0.6) 
given respectively by the blue, red and green curves. The denser the system is, the larger the energy and decay rate 
distributions are. The continuous straight lines show the timed Dicke state decay rates Tn for the three different 
system sizes. For dilute clouds (green line), the timed Dicke emission rate Fjv is centered on the distribution P{T) 
and Fat — F. When the cloud optical density increases (blue line), the timed Dicke decay rate tends to the tail of the 
distribution and Fiv ~ Fmax- 



states of the system. Some of these states iV'super) have short lifetimes Fjuper > T and are thus called super- 
radiant. Some other states jV'sub) have long lifetimes Fsub < T and are called subradiant. As the effective 
Hamiltonian is non-Hermitian, its eigenstates are not orthogonal and have common features with autoion- 
izing states or Fano resonances fTS\. Fig. |4] shows the mode energies and decay rates for three different 
system sizes. It shows that the denser the system is, the larger the energy and decay rate distributions are. 
Fig. |4]also shows the timed Dicke decay rate Fat to compare it to the decay rate distribution of the system 
modes P{V)- For dilute clouds, the timed Dicke emission rate Fjv is centered on the distribution P{T) and 
Tf4 ~ F and when the cloud optical density increases, the timed Dicke decay rate tends to the tail of the 
distribution and Vn — Tmax- 

After driving the system with the laser for a long time, which allows populating subradiant states using 
the scheme sketched in Fig. |5] we monitor the decay of the system by looking at the excited state population 
P{t) after switching off the laser. A typical decay curve of the excited state population computed from 
the numerical solution of the effective Hamiltonian approach is shown on Fig. |6] This figure shows the 
excitation probability P{i) — |/3j p as a function of time (black solid line), obtained by integrating Eq. 



www.fp-journal.org 



© 2007 WILEY- VCH Vcrkg GmbH & Co. KGaA, Wcinhcim 



12 



T. Bienaime, R. Bachelard, N. Piovella, and R. Kaiser: Cooperativity in ligiit scattering by cold atoms 



Mode couplings 

\\/Nnr) Subradiant states 

\G) -i- 



Fig. 5 (color online) Sketch of the subradiant emission couplings. 




Time (r^M 



Fig. 6 (color online) Excited state population P — J^^. \l3j\'^ decay after switching off the laser (the initial state 
corresponds to the steady state. The black dashed curve shows the single atom decay (without cooperative effects). At 
first, the fast decay corresponds to superradiance with a rate Fjv = (1 + 6o/12)r. After some time, part of the light 
remains trapped in the cloud which corresponds to subradiant emission. Parameters for the simulation: N = 2000, 
koCTR = 10, Qo ~ 0.01 r, Ao = lOr (the laser was on before during 50 F"^ to let the system reach the steady 
state). The superradiant emission diagram (blue curve) is strongly forward directed. The emission diagram of the 
subradiant states (red curve) is isotropic. The green curve shows the subradiant emission diagram averaged over height 
realizations of disorder. 



(13 1 for = 2000 atoms distributed by a Gaussian distribution with koa^ = 10. The other parameters 
are J7o — O.OIF and Aq = lOF and the laser is switched off after t — 50 F^^. The origin of time is set 
such that it corresponds to the time when the laser is switched off Under the action of the continuous laser 
excitation, the atoms reach a quasi-stationary state close to the timed Dicke state. The small subradiant 
fraction present in the atomic state after the exposition to the laser can be detected observing the excitation 
decay after the laser has been switched off The fast initial decay rate of the superradiant state is Fat = 
(1 + 5o/12)F as expected since the steady state corresponds approximately to the timed Dicke state. After 
some time, the emission rate becomes much below the single atom emission rate (black dotted line in Fig. 
|6]l. It corresponds to the subradiant emission region. At first, the subradiant decay is not exponential since 
several modes decay simultaneously. For longer times, it then ends up with a pure exponential decay, 
referred as the subradiant decay rate, when only one long-lived mode dominates |'26"28], as shown in 
the red part of the decay curve of Fig. |7] We have checked numerically the very intuitive result that the 
subradiant decay rate measured on the relaxation curves P{t) corresponds to the longest lifetime of the 
effective Hamiltonian modes | |28| . This confirms the role of cooperativity in long lived excitations in the 
cloud as investigated in ref. p4) . 

Using Eq. ( [38| ) from the previous section, we can study the emission diagram of the system. Fig. |6] 
shows the emission diagram as a function of time during the decay of the cloud. At first, the emission 
diagram of the timed Dicke state is clearly forward directed, a phenomenon reminiscent of Mie scattering. 
At longer times, subradiant modes show isotropic emission diagrams: they do not possess the symmetry 
of the laser excitation since they are not directly coupled to it. This property can be exploited in the 
experimental detection of subradiance. 
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Time (r-^) 

Fig. 7 (color online) Excited state population as a function of time for the same parameters as Fig. pl(Af = 2000, 
koo-R = 10, Qo = 0.01 r, Ao = 10 F). The laser is switched off at t — 0. The curve together with Fig. [Sjshow 
that the subradiant decay is purely exponential only after a certain amount of time (red part of the curve), when the 
mode with the longest lifetime dominates. We call this final decay rate the subradiant decay rate. In this example, 
r,,b = 8.5 10-3 r. 



In Ref. 1 26 1, we proposed to use inhomogeneous broadening schemes such as the cloud optical thick- 
ness, the cloud temperature, or the driving laser intensity as possible control parameters for subradiance. 
However, other parameters such as a far detuned speckle field, magnetic fields, or near field couplings can 
also be used. By control of subradiance, we mean two different things: controlling the population of the 
subradiant states as well as their decay rates. Exploiting these inhomogeneous broadening schemes allow 
to control and tune the dipole-dipole couplings, which is the genuine interaction leading to cooperative 
effects (superradiance, subradiance, cooperative Lamb shift). This would allow for the first observation of 
subradiant emission from a cloud of atoms in free space. 

3.3 Dipole-dipole induced suppression of excitation 

In the perturbative limit, where each dipole is driven by the external field, plus a small perturbation by the 
field scattered by all other dipoles, we can obtain an analytical solution of the excitation state population 
and the angle-resolved scattered field. 

^From Eq. ([38]l we obtain the steady-state scattered intensity in the direction {9, (p) for the timed Dicke 



state, still neglecting saturation: 



and the total scattered intensity 

r2n 



Ps = r^ I ( desmeiir,9,q)) = Po-...^. 

Jo Jo + i TV 



2 , r2 ' ^^^> 



where Iq is the incident intensity and Pq = /o/(47rfcQ). In Fig. [s] we plot the normalized excited state 
population (blue solid line on the left figure) as a function of bo and for Aq = lOOF, 

p AAl + r^ ^ AAl + r^ ^^^^ 



NPW 4A2, + r2, 4A2+r2(l + foo/12)2' 
and the normalized total scattered power (blue solid line on the left figure) 

P. r2 



(1 + 6o/12), (50) 



NPo 4A2 + r2(l + 6o/12)2 
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Fig. 8 (color online) Normalized excited state population (left) and normalized total scattered power (right) as a 
function of the optical thickness bo for Ao — 100 F. The blue solid lines correspond to Eqs. \49\ , \50\ and the red 
dashed lines to Eqs. \5l\ , ^52\ . 



where P^-^^ = ^o/ (4AQ-|-r^) is the single-atom excitation probability and we have neglected the collective 
Lamb shift, Aat w Aq. We observe that increasing the optical thickness bo the excitation population 
decreases, whereas the total scattered power has a maximum around bo ^ 24(Ao/r). The normalized 



excited state population and total scattered power resulting from the partial wave solution ( 35 i in the 
continuous-density approximation, obtained in ref. p9| for large Gaussian clouds. 



P 



4A2 



iVP(i) Aor(6o/3) 



arctan 



Aor(6o/3) 



4A2+r2 (1 + 60/6) 



and 



Ps 
NPn 



bo 



r^bp 1 + bo/12 

Y2 



4AI 



(51) 



(52) 



are also plotted in Fig|8](red dashed lines): they show a deviation from the timed Dicke approximation for 
very large resonant optical thickness 60 • 

^From these curves we can see that even for dilute clouds of cold atoms, the long range coupling 
between the dipoles leads to a cooperative modification of the excitation of the atoms and its related total 
scattered power. Note that the total scattered power initially increases with increasing number of atoms, 
despite the decrease in the normalized total population of the excited state. This can be understood by 
the fact that cooperativity leads to enhanced superradiant emission rates. For larger number of atoms, the 
suppression of the atomic excitation dominates the enhanced superradiant emission and the total scattering 
rate of the large cloud of atoms is reduced by the dipole-dipole couplings. We stress that even though the 
signature of such a suppression of fluorescence of the cloud of N atoms might bear a resemblance with a 
photon blockade regime [l^^-TSl, our model does not take into account optical nonlinearities required to 
describe such photon-photon coupling. A suppression of excitation can thus be obtained in the absence of 
nonlinear optical response. In contrast to dipole blockade effects in Rydberg states with near-field 
or Van der Waals coupling (1/r^), we are in the presence of long-range dipole-dipole couplings where 
all atoms participate in the suppression of the atomic excitation, and not only a small volume around an 
excited atom. 



4 Conclusion 

In this paper, we have presented a master equation and an effective Hamiltonian approach to describe coop- 
erative effects in clouds of cold atoms. This master equation approach, even though still restricted to single 
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excitation, allows to go beyond the effective Hamiltonian approach. In particular, we have highlighted the 
possibility of cooperative suppression of the atomic excitation, via the long range dipole-dipole couplings 
and in the absence of any non-linear photon blockade mechanism. Future work will include the possibil- 
ity of experimental observation of Dicke subradiance, long range dipole-dipole blockade and cooperative 
effects beyond linear optics. 
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